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Abstract 

An approach which unifies the Double Logarithmic Approximation at small x and the leading or¬ 
der DGLAP evolution of fragmentation functions at large x is presented. This approach reproduces 
exactly the Modified Leading Logarithm Approximation, but is more complete due to the degrees 
of freedom given to the quark sector and the inclusion of the fixed order terms. We find that data 
from the largest x values to the peak region can be better fitted than with other approaches. 


1 


The perturbative approach to Quantum Chromodynamics (QCD) is believed to solve all 
problems within its own limitations provided the correct choices of the expansion variable 
and the variable(s) to be fixed are used. However, perturbative QCD (pQCD) currently has 
the status of being a large collection of seemingly independent approaches, since a single 
unihed approach valid for all processes is not known. This is problematic when one wants to 
use a range, qualitatively speaking, of different processes to constrain the same parameters, 
for example in global hts. What is needed is a single formalism valid over the union of all 
ranges that the various pQCD approaches allow. This unification must be consistent, i.e. it 
must agree with each approach in the set, when the expansion of that approach is used, up 
to the order being considered. 

The evolution in the factorization scale of fragmentation functions (FFs) D{x,Q‘^) 
{D is a vector containing all quark FFs Dq, all antiquark FFs Dg and the gluon FF Dg) at 
large and intermediate momentum fraction x is well described |^by the leading order (LO) 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation ^ 


^ jD{x,Q^)= C ^as{Q^)P^^\z)D 


( 1 ) 


dlnQ2 

where P^^\z) are the LO splitting functions calculated from hxed order (FO) pQGD. 
We define = as/{2n), which at LO obeys as{Q‘^) = 1 /{l3o\n.{Q‘^/^ where /do = 
{11/Q)Ca — {‘^/3)TRnf is the first coefficient of the beta function and Aqcd is the asymp¬ 
totic scale parameter of QGD. For the color gauge group SU(3), the color factors appearing 
in this Letter are Cp = 3/4, Ca = 3, and Tr = 1/2; nj is the number of active quark flavors. 
On the other hand, at small x the Double Logarithmic Approximation (DLA) 


d 


;D{x,Q^) = 


dz 2 Ca 


z z 


Az 


a,(QQD(-,Q' 
z 


X 


dlnQ2 

is required, where A = 0 when D is a valence quark or non-singlet FF, while 


( 2 ) 


A = 


n 2Cf 

^ Ca 

0 1 


(3) 


when D = (P>e, Dg), where I singlet FF. The Modified Leading 



Logarithm Approximation (MLLA) improves the description here by including a part 

of the FO contribution that is known to be important at small x. With certain qualifications 
[^, the MLLA leads to a good description of all data down to the smallest x values. However, 
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what is still lacking is a single approach which can describe data from the largest to smallest 
values of x. We now construct such an approach, but leave the more detailed arguments to 
a future publication. 

As 2 ; —>• 0, the LO splitting function asP^^\z) diverges due to terms of the form ag/z. 
These double logarithms (DLs) occur at all orders in the FO splitting function, being gen¬ 
erally of the form (l/z)(as lnz)^(as In^ 2 :)'’ for r = —l,...,cx). As x decreases, Eq. ((T)) will 
therefore become a poor approximation once ln(l/a;) = The reason why Eq. (j21) 

is valid at low x is that it accounts for all double logarithms (DLs), by essentially sum¬ 
ming them up. What we want, rather, is an evolution of the form of Eq. (HD, but with the 
modihcation to the splitting function 

asP^^\z) P^^{z, ttg) + asP^°\z), ( 4 ) 

where P^^{z, Og) contains the complete contribution to the splitting function to all orders 
from the DLs, while asP^^\z) is the remaining FO contribution at LO. It is obtained by 
subtracting the LO DLs, already accounted for in P^^, from asP^^^{z) to prevent double 
counting. We will now use Eq. o to gain some understanding of P^^. For this we need to 
work in Mellin space, where the Mellin transform is dehned by 

f{uj) = [ dxx‘^f{x). (5) 

Jo 

Upon Mellin transformation, Eq. o becomes 

= 2C,a.(Q-)ADM^ ( 6 ) 

Making the replacement in Eq. in Eq. (HD and neglecting the FO term agP^^^ {z) for now, 
taking its Mellin transform 

and then substituting this into Eq. (ED gives an equation for viz. 

2(pDL)2 ^ ^pDL _ 2CAagA = 0. (8) 

We choose the solution 

Ug) = — -j- \/a;^ -|- IGCacls^ , (9) 
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since its expansion in yields at LO the result 




0 


0 


( 10 ) 


which agrees with the LO DLs from the literature 
splitting function of the 
with the results of Refs. 


. Equation contains all terms in the 
orm {as/oj){as/oj‘^Y^^, being the DLs in Mellin space, and agrees 
We now return to x space, where Eq. dH]) reads 


(^4v/^lni^ , (11) 

with Ji being the Bessel function of the hrst kind. 

To summarize our approach, we evolve the fragmentation functions according to Eq. m, 
but with the replacement of Eq. (jH), where P^^{z, a^) is given by Eq. (ITT|l . and asP^^\z) is 
given by agP^^Yz) after the terms proportional to Ug/z have been subtracted. 

Before we outline the phenomenological investigation of our approach, we note that it is 
more complete than the MLLA, which can be shown as follows. With asP^^\z) accounted 
for, Eq. Q is modified to 

+ ^o; + o,s{Q^)P^ QY-, 

up to terms which are being neglected in this Letter and which are neglected in the MLLA. If 
we approximate asP*'°\a;) by its single logarithms (SLs), dehned at LO to be the coefficients 
of a;°, equal to those in asP^'^\uj), 




0 


- 30 ? 


iTnUf -^-Ca- ITrUj 
then if we apply the approximate result that follows from the DLA at large Q, 

Dgrq = 


(13) 




(14) 


(e.g. this can be derived from Eq. (CH)), the gluon component of Eq. m becomes precisely 
the MLLA differential equation. Therefore we conclude that, since we do not use these two 
approximations, our approach is more complete and accurate than the MLLA. 
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We now test our approach by comparing its effects on fits of quark and gluon FFs to data 
to the standard FO DGLAP evolution. We use normalized differential cross section data 


for light charged hadron production in the process e'''e —>■ ( 7 , Z) —h + X, where h is the 


veu nauion ana is anyiimiK eis e, iiom ai ys = 

Q, MARK II |nl and TPC Q at 29 GeVJ^OPAZ at 


GeV 

DELPHI Q, L3 


Q, OPAL 


GeV j23[. We place a small x cut on our data of 


GeV 


a 


observed hadron and V is anything else, from TASSO at y/s = 14, 35, 44 GeV |9j] and 22 

,^LEPH a> 

19l[ an d OPAL 
22 1 and 202 


and MARK II [18| at 91 GeV, ALEPH 


at 133 GeV, DELPHI at 161 GeV j2l| and OPAL at 172, 183, 189 GeV 


^ = ln(l/x) < 


(15) 


where M is a mass scale of 0(1) GeV. We ht the gluon g[x, Ql), as well as the quark FFs 


fuc{x, Ql) = ^ (m(x, Ql) + c{x, QD) , 
fdsb{x, Ql) = ^ (d(x, Ql) + s{x, Ql) + h{x, Ql)) , 


(16) 


where Qo = 14 GeV. Since the hadron charge is summed over, we set Dg = Dg. For each of 
these three FFs, we choose the parameterization 


/(x, Qo) = ^ exp(—cln^ x)x"(l — x)^. 


(17) 


which at small x is a Gaussian in ^ for c > 0 with centre positive in ^ for a < 0 as is 
found to be the case, while it reproduces the standard parameterization (i.e. that without 
the exp(—cln^x) factor) used in global hts at intermediate and large x. We use Eq. fimi to 
motivate the simplihcation 


^uc ^dsb 1 

(18) 

^uc ^dsh 

to our parameterization. We also ht Aqcd, giving 9 free parameters. Since we only use data 
for which y^ > rrib, where rrib ~ 5 GeV is the mass of the bottom quark, and since Qo > mb, 
we will take uj = 5 in all our calculations. While the precise choice for Uf does not matter 
in the DLA, calculations in the FO approach depend strongly on it. 

If the (1 — x)^ factors were absent, Eq. (HI would dictate that 

Nuc = N^sb = ^Ng. (19) 

La 
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TABLE I: Parameter values for the FFs at Qq = 14 GeV parameterized as in Eq. (HZD from a 
fit to all data listed in the text using DGLAP evolution in the FO approach to LO. Aqcd = 388 


MeV. 


Parameter 

FF 

N 


a 

c 

g 

0.22 

-0.43 

-2.38 

0.25 

(u+c)/2 

0.49 

2.30 

[-2.38] 

[0.25] 

(dH-sH-b)/3 

0.37 

1.49 

[-2.38] 

[0.25] 





FIG. 1: Fit to data as described in Tabled Some of the data sets used for the fit are shown, 
together with their theoretical predictions from the results of the fit. Data to the right of the 
horizontal dotted lines were not used. Each curve is shifted up by 0.8 for clarity. 

However, the (1 —factors are important at large x, and Eq. (j14jl is only an approximation 
at small x. Thns it will be interesting to see the deviations from Eq. (HI after htting. 

We hrst perform a £t to all data sets listed above nsing standard LO DGLAP evolntion. 
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i.e. Eq. ([T]) without the replacement in Eq. (jH). We ht to those data for which Eq. (ITH|l is 
obeyed with M = 0.5 GeV. This gives a total of 425 data points out of the available 492. 
We obtain Xdf — 3.0, and the results are shown in Fig.[T]and TableHl The result for Aqcd is 
quite consistent with that of other analyses, at least within the theoretical error of a factor 
of 0(1). It is clear that FO DGLAP evolution fails in the description of the peak region and 
shows a different trend outside the ht range. The exp(—cln^a:) factor does at least allow 
for the ht range to be extended to x values below that of a; = 0.1, the lower limit of most 
global hts, to around x = 0.05 = 3) for data at the larger a/s values. Note that j3g is 

negative, while kinematics require it to be positive. However, this clearly does not make any 
noticeable diherence to the cross section. 


TABLE II: Parameter values for the FFs at Qo = 14 GeV parameterized as in Eq. m from a fit 
to all data listed in the text using DGLAP evolution in our approach. Aqcd = 801 MeV. 


Parameter 

FF 

N 

d 

a 

c 

g 

1.60 

5.01 

-2.63 

0.35 

(u-|-c)/2 

0.39 

1.46 

[-2.63] 

[0.35] 

(dH-sH-b)/3 

0.34 

1.49 

[-2.63] 

[0.35] 


Now we perform the same ht again, but using our approach, i.e. Eq. m with the replace¬ 
ment in Eq. 0, for the evolution. The results are shown in Table Inland Fig. |21 We obtain 
Xbf — 2.1, a signihcant improvement to the ht above with FO DGLAP evolution. This 
should also be compared to the ht to the same data in Ref. |^, where DL resummation 
was used within the MLLA but with neither FO terms nor quark freedom (i.e. Eq. (HI was 
imposed over the whole x range) and Xdf = ^.0 was obtained. The data around the peak 
is now much better described. The energy dependence is well reproduced up to the largest 
a/s value, a/s = 202 GeV. We conclude that, relative to the MLLA, the FO contributions 
in the evolution, together with freedom from the constraint of Eq. dH, make a signihcant 
improvement to the description of the data for ^ from zero to just beyond the peak. How¬ 
ever, Aqcd is rather large, even within the theoretical errors. Ng is too large by a factor 
of about 2 relative to its prediction provided by Eq. da. However, note that Ng is weakly 
constrained since the gluon FF couples to the data only through the evolution, requiring 
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FIG. 2: Fit to data as described in Table ITIl 


e.g. gluon data to be properly constrained. These problems are related to the worsening 
description of the data on moving beyond the peak, since hts in which the cuts were moved 
to larger ^ values gave an increase in Aqcd and Ng, as well as Xdf- Figure El is repeated 
in Fig. El to show more clearly the good quality of the £t at intermediate and large x. A 
couple of points at x = 0.9 are not well described, although the data here are scarce and 
have larger errors. 

In conclusion, we have proposed a single unihed scheme which can describe a larger range 
in X than either FO DGLAP evolution or the DLA. Further improvement in the small x 
region can be expected from the inclusion of resummed SLs. Alternatively, improvement 
may be achieved by suppressi^ the higher moments’s evolutions, since these are unstable 
yet formally of higher order |2|, and the suppression of these effects provided by the FO 
contribution is unlikely to be sufficient. Our scheme allows a determination of quark and 
gluon FFs over a wider range of data than previously achieved, and should be incorporated 


into global fits of FFs such as that in Ref. 
limited. 


since the current range of 0.1 < x < 1 is very 
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■ OPAL 91 
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• TPC 29 

o TASSO 14 



FIG. 3; As in Fig. [2 but with the cross section on a logarithmic scale versus x. Each curve, apart 
from the lowest one, has been rescaled relative to the one immediately below it by a factor of 5 for 
clarity. 
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